LIMIT THEOREMS FOR COMPETITIVE DENSITY-DEPENDENT 

POPULATION PROCESSES 
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1. Notation 



We write No for the set of non-negative integers. 

We use Landau asymptotic notation, where all asymptotics are with respect to N, defined 

below. We write f(N) = 0(g(N)) and f(N) = o(g(N)) if 

f(N) 



and 



lim sup 

7V->-oo 



lim 

JV->oo 



f(N) 



g(N) 



< oo 







respectively. We also use Hardy notation: f(N) <^ g(N) if f(N) = o(g(N)). 
Throughout, I will use d{ to indicate the partial derivative with respect to the i cooordinate 
and D to denote the total derivative operator: if F : IR n — > W 71 , 

(DF) = (djFfo 

Rf = {x g R K : xt > 0, i = 1, . . . , K}. 



If x, y E M , we write x < y if Xi < yi for all i, x < y if x < y and X{ < yi for at least one 
i, and x < y if x% <y% for all i. 

Given a Polish space fi, Dq[0, oo) denotes the space of £1- valued cadlag functions endowed 
with the Skorohod topology. 

We use X — > X to denote convergence in distribution for a sequence {X (£)} of cadlag 
stochastic processes, i.e. 

E[f(X N )]^E[f(X)] 
for all / : B^fO, cxd) — >■ R continuous in the compact uniform topology. 



2. Introduction 



3. Competitive density-dependent population processes 



To begin, I will introduce the object of study, intended to encompass a variety of models considered 
in population dynamics, population genetics, and community ecology. Throughout, I will consider 
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a population consisting of a finite number of discrete individuals. Each individual has a type, which 
is indexed by an integer in the set {!,... , K}. All individuals of a given type are identical. 



Denote the number of individuals of type i G {1, . . . , K} at time t by Xf (t) and let 
I assume that the number of individuals of each type changes when there occurs 



'1) a reproduction event in which some individual of type i produces a clutch consisting of 



def 

|n| = n\ + • • • + tik offspring, of which n,- are of type j, or 



(2) the death of an individual of type i. 

Thus ~K N (t) is a Markov chain on Ng , which may be represented as 

X^(t)=X» + f; £ nB^-^Dfit) 

j=l neN^ »=1 

where B^ n (t) and D^(t) are counting processes with intensities 

#n(^K(<) and *f(^)*f(*), 
respectively, i.e. 



(i) <„w = <„(*) - j Pi!n (^ J*r (s) ds 

and 



'0+ 

are martingales, with quadratic covariations 



f(^)x^s)ds 



Jo+ 
covariations 

[^] t = <W and [Df] t = Df(t) 



In this work, we focus on identifying and providing rigourous proofs of limiting processes. Elsewhere, 
we discuss in detail the biological implications of applications of this approach to population genetics 



3.1. Assumptions. I will assume that for all compact sets % C R , 

(2) 2, l n l SU P A,n( x ) < °° an d SUp 2, l n l su Pft,n( x ) 



and that 



nSN£ neN/5 



(3) lim N sup /3f ne . (x) - /3 in (x) = and lim N sup 6? (x) - J< (x) = 



uniformly on compact sets, while for all i,j and all compact % C R , 

oo 

(4) ^2 supn/3 in (x) < oo 

and 



n=l 



xG3C 



(5) 



sup 

xeuc 



E 

{neN^ : n^|n| ei } 



n^ n (x) = 0(±) 



Thus, with probability tending to 1 as N — > oo, all offspring are of the same type as the parent, 
i.e. mutation is rare. 



While these summability assumptions are made for technical reasons, they are eminently plausible 
from a biological perspective, as they simply requires that the mean and variance in the number 
of offspring produced in any single reproductive event are both finite. Both are readily satisfied by 
assuming some fixed maximal clutch size. 



3.2. Mutation. In this section, I will introduce new notation that allows the birth-death-mutation 
process to more closely resemble the forms considered in the Wright-Fisher diffusion. Let 

^(x) = V |n|/C(x) 



n6N^ 



and 



N I \ dcf I 



<(x) 



/?f(x) 



E n Aw- 



{nGN ( f:n^|n| ei } 



f3j (x) is thus the expected total reproductive output of an individual of type i per unit time in 



environment x, while /U^ (x) is the fraction of the expected number of offspring which are of type 
j. By assumption, 

4(x) = 0(^). 



4. Law of Large Numbers 



Let 

and define F : R^ -)• R^" by 

Let tpt denote the flow of 

(6) ^^^^^ 

and let 



Pi(x) = ^nft„( x ) 



n=\ 



clef 



Fi(x) = (A(x)-<Ji(x))xi 

Y(i) = F(Y(i)) 

Y"(t) = £X"(t). 



We then have the following Law of Large Numbers for Y (£), which will be proven below in Section 



Proposition 1. Let x € R+ , and fixO<r<s<l and e > so that 

def ■ 



X X)£ - {y : lly - <^ x ll < e f° r some * > 0} C int(R^). 



Then there exists a constant B^ £ such that 



sup llY^^-^xll >N~ r 



1-s 
3 x,e 



0<-S^-lnAf 

— iBi 



N~ r 
> < 



In particular, for all t > 0, 



lim Y JV (t) - iptx \\= P x -a.s. 



Based upon this, I will say that the Markov process X Ar (t) is competitive if the dynamical system 
^j is competitive [5], i.e. /3j(x) and 5j(x) are C 1 and 

dj (A(x) - <^(x)) <0 

for all j ^ i. 



4.1. Competitive Dynamical Systems. Subsequent to Proposition [H I will always assume 
X (i) is competitive, and further that the dynamical system ([6]): 



(1) is dissipative: there is a compact set X that uniformly attracts each compact set of initial 
values, 

(2) is irreducible: the matrix (dj (/3j(x) — <5j(x))) is irreducible for all x € intR+, 

(3) has a source at the origin, and 

(4) for i y^ j, dj (/3j(x) — <5j(x)) < at every equilibrium in R^ \ 0. 

Again, each of these assumptions has a biological interpretation: the first requires that the pop- 
ulation remain finite, the second that all types interact with all other types, and the third that 
populations will grow when started from small initial densities. The last is required 



5. Linear Birth-Death Process Approximation 

f 

6. Quasi-Neutrality and Weakly Selected Quasi-Neutrality 



Henceforth, we will say that the process X Ar (t) is quasi-neutral if there exist C 2 functions 7 : R^ 
Rf and R : R+ -> R such that 

(7) FiW^-yiWRW. 



and 



(1) 7(0) >0, 

(2) R(0) = 1, and {x G Rf : i2(x) > 0} is compact, and 

(3) (DR)(x) < for all xGlf. 



Under these assumptions, 

Q = {x : R(x) = 0}, 

is an attracting, compact, co-dimension one C 1 -submanifold of equilibria for the dynamical system 
([U]), as can be seen by considering the Lyapunov function V(x) = i? 2 (x). 



6.1. Geometry of fL f2 is diffeomorphic to the standard simplex 

A*" 1 d =/ |x € 
via the radial projection map 



i=i ) 



P(x) = — R x. 

2^i=i x i 



diPk = 7^7^ ($ik — Pk) 



This has partial derivatives 

1 

E/=i xi 

and 

didjp k = — -J (§ik + 5 jk - 2p k ) 

(l£i *i) 

and is thus differentiable when restricted to Q, where the denominator is non- vanishing. 

To see that this is invertible, consider f x (t) = R(tx) and F(x, i) = (x, f x (t)). Then, 

det((Z>F)(x,i)) = (d t f x )(t) = (DR)(tx) • x, 

which, by Assumption ([2]) is strictly positive. Thus F~ l exists and is differentiable; in particular, 
/~ 1 (0) exists, is unique, and is differentiable in x, so x 4 / x r 1 (0)x : A K ~ X — > S7 is a differentiable 
injection. 

We will henceforth write 

a naming convention whose origin will be discussed in greater detail below. 

5 



6.2. Local Dynamics Near Cl. Recalling ((7J), we see that 

(djFi)(x) = (d j ji)(x.)R(x)x i +ji(-x)(djR)(x)xi + 7i(x)i2(x)<%. 

In particular, writing 3~i(x) = 7j(x)xj, if x* G fi, 

(DF)(x*) = F(x*) ® (Di2)(x*), 
Thus, (DF)(x*) has two distinct eigenvalues, and A(x*), where 

(8) A(x) d ^ f (Di?)(x).?(x) 

The former has corresponding eigenspace (Di?)(x*)- L = T x *Cl, whilst the latter corresponds to the 
single eigenvector 5F(x*). 



6.3. The Projection Map and a Time Change. Let tpt denote the flow of ©, and let 

7r(x) = lim ipt*-, 

t— >oo 

POO 



T x (t) d = J°° R{xl, 8 y)ds, 



it 
and 

r(x) d ^ f T x (0). 
Prom the definition of T x , we see that 

7^ x {t-u)=7 x (t), 
whence r( , 0jx) = T x (t), a relation we shall need later. 



7. Results 



7.1. Quasi-Neutral Processes. Using these definitions, we show that Y^ will approach O-limit 
set with probability tending to 1 as N — > oo: 

Proposition 2. Let e and B Xj£ be as in Proposition^ Fix r < —5 — . Then for all r' < r < 

i+ /3_(l-a) 

s < 1, and N sufficiently large, 



r[Y N (j^lnJV 



> N--} < * 



£ 



This establishes that the for N sufficiently large, Y (i) will come arbitrarily close to Cl. In fact, 
the process, having arrived at Cl, remains there. To do show this, I consider the process on a longer 
time-scale: let 

Z JV (i) def JL X "(ft). 

Then, 



Theorem 1. Let r' be as above, and fix < 5 < r' . Suppose r(Z (0)) < N , and let 



T-Sl 



rs = M{t > : T(Z N (t)) > N~ d }. 



Then, as N — > oo ; for any fixed T > 0, 



t(Z n (TAt s ))-^0. 



Recalling that z = (f) T ? z \7r(z), the continuous mapping theorem gives us 



Corollary 1. Under the assumptions above, 



Nu\\ B > 



Z N (t) - ir(Z N (t)) ^U 



We now prepare for our main result by introducing some additional notation. Let 



(9) 



<Tj(x) = lim N 



/3f(x)-A(x) -(^(x)-ft(x)) 



i.e. <Tj(x) is the 0(]y) component of the net reproductive rate (that the limit exists is guaranteed 
by Assumptions (J2J) — ((5]), let 



(10) 



9y(x)^ lim iV4(x) 



be the rescaled rate of mutation, and let 



A(x) d ^^n 2 ft in (x) 



n=0 



With these, we have 



Theorem 2. Suppose that as N — >■ oo, Z Ar (0) — > Z(0) € f2, and that the Martingale problem for 

K / K A 

(11) (Lf)(n) = £ -J>i(ir)ft(ir)7r« + EM*)&(*)^ 

i=l V 7 = 1 7 = 1 



X A 



+ ^r^r E E ^ « w») - (W*)) ^ 

^ 7rj i=i fc=i 



A" 



* I **(*) - I^E^^w^w^ 



+ ^y [2 J>ili) ((a,-iJ)(ir)(^(7r) + &(»)) - (^)(7r)(A(7r) + &(*))) 7i ( w ) ffi 



i=i 



A" A 



+ is w E Ewww^x*)^ ^TT ^» + & w)w 



i=l fc=l 

A 



7j(t) 



7»(t) 



+ 



bww ^^y j (4-w+^-wk, 



A( 

AAA 



EEX>^)Vh- ( (^(I) ) (ftw+ftW)^7 fc wwi(^)^ ) ) )(^/)(7r, 



A A 



+ EEM (aw+a(t))^ 



7T; 



A 



^(9 fcj R)(7T) 



A 2 (tt) ^ 
1=1 J=l \ v I fc=l 

+ 7j(w)7fc(7r)(0ii2) (A(w) + ft(ir)) - 7 i (7r) 7i (7r)(5 fej R) (&(ir) + &(*)) 



7i(w)7fc(w)(a,-R)(^(7r) + ^(ir)) 

vfc \{didjf){iv) 



is well posed. Then II (t) = ir (Z N (t)) converges weakly to a diffusion process H(t) with this 
generator. 



In general, well-posedness for degenerate diffusions remains an open-problem, and uniqueness is not 
guaranteed when the coefficients fail to be Lipschitz; in section 19.5.51 1 discuss sufficient conditions 
to ensure uniqueness. 



Note that the assumption that Z (0) — > Z(0) £ Q, is necessary to obtain weak convergence in the 
results above; Propositions Q] and [2] tell us that in general, 



lim Z"(0+) = 7r(Z(0))#Z(0). 

N— >oo 

Thus, the limiting process fails to satisfy the Feller property. We may still obtain a global conver- 
gence result, by considering a variant of the original process: 



Corollary 2. Let 



rN V 



Z N (t) = Z N (t) - ^ 4 Z» + 7r(Z"(0)). 



Then Z — > Z, a diffusion process with generator (jllj) as above. 



In keeping with the convention in population genetics of considering relative frequencies, rather 
than absolute numbers of types, we consider the process P(t) = p(Jl(t)), where Applying Ito's 
formula to p(II(t)) yields the following: 

Corollary 3 (Relative Frequency process). 

(12) L RF f( P ) = W -xx-(p)#(p)tt +E fi(p)^(p)» 

»=i V j=i i=i 

+ (7f(p) A 7 p ; (p))K £ £ %(p) ((ww - ( w(p)) to 

+Pi ( K(p) - ^ p (p)) - (7f(p ^J (p)) £(w(pK(pM- 



(7f( P )-r(p))^^^ fl mw » Jfl aw(p 



(a^(p)) 3 t^tT"'"' "" ~ j V(W(p 



j=i k=i i=i 



Where we adopt the notation f p (p) = f(p 1 (p))- 



Lastly, we observe that when all types are identical to order 0(-^), i.e. 

7i (x) = 7 3 -(x) for all 1 < i,j < K and all x € M+, 

we obtain the following process, which naturally generalises the Wright-Fisher diffusion and the 
diffusion approximation to Gillespie's fecundity- variance model OH]: 

Corollary 4 (Weak Selection). 
(13) L WS f(p) 



8. Important Special Cases 



8.1. The Gause-Lotka-Volterra Model. 



8.2. The Double Monod Model. 

/ 

8.3. Results for K = 2. 



9. Proofs 



9.1. Proof of Proposition [Q Define F^ : Rf -> R^ by 



Thus, 



K K 



X N (t) = X» + r iVF^(^)) ds + £ £ nBf n (t) - ^e^t) 

,/0+ ' »=1 nGN ( f *=1 

and 



Let 



i=i neN K ■ i=i 



A,^ = sup ||y|| . 
y60C x , e 



Next, recall that F is C 2 and thus locally Lipschitz, so there exists a constant -B Xi£ such that 

||F(yi)-F(y 2 )|| < B^ E || yi - y 2 || 
for all yi,y 2 € 3C X , £ - 

Now, define a stopping time 

rf = inf{£ > : || Y N (t) - ^x|| > e} 
and consider the stopped process Y N (t A r^ ). 
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From ??, we have 



tAri' 



Y N (t A rf ) - Y{t A rf , x) = Y N (0) - x + / F(Y JV (a)) - F(V> s x) ds 

+ r TiV F^(Y^( S ))-F(Y^( S ))d S + l (f; £ n^^Ar^-^e.A^Arf)] , 
0+ \i=l neN( f i=l / 



while applying the Cauchy-Schwartz inequality gives 



\Y N (t A rf ) - Y(t A rf,x)|| 2 < 4 II Y w (0) - xl 



+ 



tAr£ 



F(Y lv (s))-F(i; sX )d S 



0+ 



+ 



tAr£ 



Ni 



F n (Y n (s))-F(Y n (s))ds 



0+ 



+ 



N 2 



K 



K 



E E ^A^-Ee^Ar, 



e i 



*=1 n£N ( f 



j=l 



Applying Jensen's inequality gives 



and 



sup 



sup 



tAr£ 



tAr£ 



F(Y JV (s))-F(^x)ds 



0+ 



< Bi £ T 



TAt£ 



Y n (s)) - ip s x\\ ds 



0+ 



N / 



F n {Y lv (s))-F{Y lv {s))ds 



<T 



j-TAri 
J0+ 



\F N (Y N (s)) - F(Y N (s))\\ 2 ds, 



while 



F t N (y) - Fi(y) = £ n (/^(y) - /3 m (y)) 



yi{s) 



n=l 



K 



+ E E 4(y)»;(') - (*f(y) - *(y)) ^( s )> 

i=l {neN^:n^|n| ei } 

so by assumption, we have ||F^(y) — F(y)|| = O(^) uniformly on 3C X)£ . 



Lastly, Doob's inequality gives 


E 


sup 

t<T 


E E ^J> A rf ) - L 

3=1 n£N ( f 

<N E 



rf; 



<4E 



K 



EE*-^f 



i=i neN ( f 



(T A rf ) 



A' 



E E n l ^{Y N {s))Y l N {s) + 5f{Y N {s))Y l N {s) 



i=l nGNjf 



(is, 



11 



and again, assumptions made on the intensities insure that 
uniformly bounded on X* £ . 



^K 



cNi 



Ej=i EneN^ n iPj,n(y)yi + 5 i (y)v 



IS 



Combining these, we have 



E 



sup llY^t A rf) - Y(i A rf )|| 2 



t<T 



<4E 



Y JV (0)-Y(0) 



+ %r + 4i ^ T 



TAri' 



E 



for a constant C X)£ . Applying Gronwall's inequality, we have 



(14) E 



t<T 



AMI 2 



sup Y JV (tArf)-Y(tArf) 



< 4E 



Y JV (0)-Y(0) 



Y"(s)-Y(s) 



C x , e T \ 4B 2 T 

+ N l ' 



ds 



from which, taking T = } B z In N, we obtain 



sup E x ||Y Ar (tAri v )-Y(tAr £ Ar )|| 2 < Cx ' e(1 "° 



t<^^-lnJV 



4B» 



(In N)N~ S <N- 



for iV sufficiently large. 



Lastly, we observe 



{tArf <*}= P x { || Y^ (t A rf )- Y(i A rf)f> £ } 



< 



E x Y^tArf )-Y(tAr f 



-AMI 2 



< 



AT" 



from which our result follows. 



9.2. The Projection Map and its Derivatives. We begin by considering a new dynamical 
system, with trajectories identical to ([6]), but traversed backwards in time, away from Q. 



Lemma 1. Let (fit be the flow of 

(15) f ^^ 

then 

(16) 



y = -^(y) 



^tx = (/> Tx(4 )7r(x). 



Proof. Recall that 
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First, assume that x is a rest point for ([6]); then, for all t, iptx = x and R(iptx) = 0. Thus, T x (t) = 
and 

V> t x = x = tt(x) = Tx(i) 7r(x) 

as desired. 



If x is not a rest point, then R(iptx) is strictly positive or strictly negative for all t£l, and thus 
the inverse function T" 1 is everywhere defined and has derivative 



■k'M - 



dt" x v " y R(ip 7 -i , t ,x) 

Thus 



while 



limi/vr- i/*iX = lim V^x = 7r(x). 



The result then follows from uniqueness of solutions. □ 



Henceforth, I will use (fit in favour of ipt when considering the trajectories of ([6]). In particular, 
taking t = in (I16D gives the following essential identity: 

(17) x = </> T(x) 7r(x) and tt(x) = 0_ r(x) x, 

As a first application, this may be used to find the derivatives of 7r(x) and t(x). Differentiating 
(fTTD yields 

(Dtt)(x) = e -/o7 (x) ( D3r )(^ x ) rfs + £F(tt(x))(Dt)(x) 

where we have exploited the identities [6] 

^Ux = -2-(&x) and (D0 t )(x)=e-C( D3r )(^ x )^. 



Differentiating R(7v(x)) = 0, we have 

= (D J R)(7r(x))(D7r)(x) 
= (Di?)(7r(x)) (e-^+ r(x) ( D3r )(^ x ) ds + 5F( 7r (x))(Dr)(x)) 
whence 
(18) (Dr)(x) = - I J^(Di?)( 7 r(x)) e -/o7 (x) (^)(^-)^, 

Substituting this into the expression for (D7r)(x), above, we have 

(Dtt)(x) = (l- x ^-i-^y (7r (x))(Di?)(7r(x))) e -/o7 (x) (^)(^-)^ 
= fl _ J_(DF)( T (x))) e -/ T w (D?)(^)^ 

V a(7t(x)) ; 

13 



(19) 



To find the second derivatives, recall that the matrix exponential has Frechet derivative [T] 

A+hB _ C A r l 



so that 

-too 



l im £ —=\ e uA Be^-^ A du, 

HO h J 0+ 

f e u /o+ T(x) ( D3r )(^ x ) ds (-(D^)(7r(x))(a,r)(x)+ f "* d t [(DJ)(^x)] ds) e (1 ~ u) ^ ^^ ds du. 

Some matrix manipulations give 

d % (A(tt(x))) = (y(7r(x)) i (D 2 J R)( 7 r(x)) + (DJl)(ir(x))(DS)( W (x))) (^7r)(x) 



K K 



J] ^ ((5 Hj R)(7r(x))J fe (7T(x)) + (0 fc i2)(7T(x))(5 / J fc )(7T(x))) ($7T,)(x) 



fc=l Z=l 



( ^ DT)(X) = (- ^ A (l ff (x)f (DiZ)(,r(x)) + ARx)) ^^^W) e/o7 (x) (D^)(^x), s 



A(7T(XJJ V 



(^D7r)(x) = $ ^/o7 (x) ( D3r )(^ x ) ds ) - (D^)(7r(x))(a i 7r)(x) <g> (Dr)(x) - ^(x)) 8) (5,Dt)(x) 
/a,(A(7r(x))) 



A 2 (vr(x)) 



(DF)(tt(x)) 



x e 



+ ^^y ((Dy)(7r(x))(^7r)(x) ® (Di?)(7r(x) + <F(tt(x)) ® (D 2 «)(7r(x))(d,7r)(x)) 

The various partial derivatives can then be obtained via the identities (c?jt)(x) = (Dr)(x)ej, 
(diTTk) = e^ • (D7r)(x)efc, et cetera. 

To end this section, we develop several other identities that will be useful: first, note that 

^U(&x) = -(Di2)(&x) • ^(^x) = -A(^x), 
ere 

so that 



(20) fl(^x) = i2(^r x(t) 7r(x)) - i?(7r(x)) = - / A(0 s tt(x)) , 

J0+ 



, ids, 
'o+ 

and, in particular, 

(21) fl(x) = - / A(</> s 7r(x))ds. 

J0+ 
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Secondly, 

j t n<t> t x) = -(dj)(^)^x), 

whence 

(22) y(&x) = e" /o+( Dy )(^ x ) <^( x ). 



9.3. Proof of Propostion [2j By the defintiion of T x (i) and ([20]) . we have 

<9 t T x (t) = -fl(^x) = / A(^tt(x)) ds < A x , £ T x (i). 

Jo+ 

where 

A X>£ = sup A(y) < 0. 

Applying Gronwall's Inequality gives 

T x (t) < e A -*T x (0), 
or, equivalently, r(^x) < e x > e V(x). Taking £ = . 1 ^ s lniV, we have 

X,£ 

r(^x) < A^ Ax ' Eil ^r(x). 
Lastly r(x) is C 2 and thus locally Lipschitz, so 

\r(Y N (t)) - r(^x)| < L \\Y N (t) - &x|| . 

Thus, provided r < tt — i+f 2 ^ anc ^ r ' < r < s < 1> then — r < A Xj£ ^^-, and, if N is sufficiently 
large, 

r(Y%JgL In JV))| >*"'■' 

only if ||"Y" JV (t) — V't x || > iV _r . The result then follows from Proposition [TJ 



9.4. Proof of Theorem [TJ. Our strategy will be to show that r 2 (Z (i A r^)) converges weakly to 
0. Then 

, {t A „ < t , _ , {Az> A „„ > , } < itW^rii - o 

o 

as N — > oo. This allows us to conclude that r(Z (£)) — 5- 0. We begin with a few preliminaries 
that are essential to our argument. 

First, we recall that we have assumed that A(x*) < for all points x* E Q. Since SI is compact, we 
can fix 5 > such that A(x) < for all x € 0$, 

ft/^ (J «mo). 

tS[-<5,<5] 

By virtue of the joint continuity of ^x in £ and x and the compactness of $7, fi^ is a compact set. 

Let 

t s = inf {t > : |r(Z"(a))| > 5} 
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/ 

Then, (fTTj) gives 

Z 7V (t A T 5 ) = cP T{Z N {tATs)) 7v(Z N (t A t s )) € fi* 

i.e. Z Ar (t A Tg) is compactly contained for all t > and all iV. 



Next, as before, we have the following integral equation for Z N (t): 

(23) Z N (t) = Z N (0) + f NF N (Z N (s-)) ds + M N (t), 



where 

M N (t) d ^ 1 ^ nB^(Nt) - J2*D?m 

is a square integrable martingale with component-wise quadratic variations 






while for i 7^ j, [M- ■ , M- ]t = 0. and corresponding Meyer process 

(24) 



[Mf) t = f £ n 4 /3f n (Z^( S -))Zi v ( S -) + ^ ( 5f(Z^( S -))Zi v ( S -) ( i S . 



In particular, we recall that by our assumption, sup^- eN SneN K n /^i^i(^ Ar ( s )) an d su PiVeNo <^i (Z (s)) 
are uniformly bounded on fij. Thus, for t < t$, there exists a constant C £ such that 

(25) ^^^^^ (M?) (t) < C £ t. 

We now turn to the core of our argument. Since Z N (t) is a quadratic pure jump process for all N, 
Ito's Formula takes the following simple form when applied to r [10] 

(26) r(Z N (t)) = r(Z») + J2 f (d l r)(Z N (s-)) dZf{s) 

+ J2 \t(Z n (s)) -r(Z N (s-)) -J2(diT)(Z N (s-))AZ?(s)\ 

0<s<t *■ i=l ' 

Formally expanding the he right hand side of (|26p in powers of N, and recalling that F (x) — F(x) = 
0(jr) uniformly on on compacts, we see that the highest order term is 

N I (Dt)(Z n (s-)) -F(Z N (s))ds 

= ~N f x( , 7 ^ M ,. (D^)( 7 r(Z^( g )))e/o7 (Z " W) (^)^^W)^(Z^( S ))^(Z^( g ))& 

J0+ A(7T(Z "(S))) 
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which, using (j22|) simplifies to 



-N I , , ^ N{s))) (VR)(7v(Z N (s)))?(7v(Z N (s)))R(Z N (s))ds = -N J* R(Z N (s))ds 



t 

o+A(7r;Z' 



Lastly, from (|2Tj) . this is 

N / / \(<j) u ir(Z N (s)))duds, 

For t < Ts and N sufficiently large, there exists A + > such that 

A(0 u7 r(Z Af ( S ))) < -A+ < 0, 
so 

nt ft r (Z N (s)) 

N X + T(Z N (s))ds<-N / \(0 u Tr(Z N (s)))duds. 

We proceed by imitating the proof of linearised stability [6j, where —N\ + acts in the role of an 
upper bound to the eigenvalues of the Jacobian, To this end, consider 

e 2NX+t r(Z N {t)) 2 . 

Again using Ito's formula to this function, we have, after some simplification, 

(27) r(Z N (t)) 2 = e- 2NX+t T(Z N (0)) 2 

ft ( ft fT(Z N ( S )) \ 

+ 2Ne 2NX+( - s - t) [X + t(Z n (s-))+ / \(<p u Tv(Z N (s)))duds )r(Z N {s-))ds 

Jo+ \ Jo+ Jo+ J 

+ I* 2Ne 2NX+ ^T(Z N (s-))(DT)(Z N (s))-(F N (Z N (s))-F(Z N (s))) ds 

+ f 2e 2NX+ ^T(Z N (s-))(d l T)(Z N (s-))dM N (s) + e N (t), 



10+ 

where 



(28) e N (t)= J2 e 2NX+ ^- t ^r(Z N (s)) 2 -r(Z N ( S -)) 2 

0<s<t *• 

K 

- ^2r(Z iV ( S -))^r)(Z JV ( S -))AZf( S ) \. 



i=l 



By our choice of A + the second term on the right is always non-positive, so 

supr(Z N (t At 5 )) 2 < e- 2NX+T ^r(Z N (0)) 2 
t<T 



+ / 2Ne 2NX+ ^-^T(Z N (s-))(Dr)(Z N (s)) ■ (F N (Z N (s)) - F(Z N (s))) ds 

Jo+ 

+ / 2e 2NX+ ^~ T ^r{Z N {s-))(d i T){Z N (s-)) dM N (s) + e N (t). 

Jo+ 



By assumption the first converges in distribution to 0. Theorem 4.1 in [2] tells us to show weak 
convergence to 0, it suffices to show that the remaining terms converge to in probability. 
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We have assumed that F (x) — F(x) = 0(-/\f) uniformly on on compacts, so the third term is 
bounded above by a constant multiple of 



,2NX+{ S -t) ds 



' 'l_ e -2iVA+* 



0+ 



2iVA+ 



The fourth term is a square-integrable martingale with Meyer process 

N 



f (2e 2NX+ ^T(Z N (s-))(d t T)(Z N ( S -))) 2 d(M N )(s). 



We recall Lenglart's inequality [8j: for < p < 2, there exists a constant C p such that 



E 



N( n \\P 



sup \M iy {s 

t<TAr 6 



<C P E 



M N ) 2 (TAt s ) 



In particular, when p = 1, the above, combined with Jensen's inequality gives 



E 



sup \M N (s)\ 

t<TAr s 



<C p E[(M n )(TAt 5 )]*. 



Applying this, and remembering that r, d%T are differentiable, and thus uniformly bounded on 
compacts, we see that the third integral is bounded by a constant multiple of 



E 



TAts 

AT I e 2NX+ ^- TAT6 U(M N )(s) 

o+ 



which, by ([25]) . is bounded above by a constant multiple of 

rT 



J2N\+(s-TAt s ) ds 



0+ 



Lastly, Taylor's theorem gives 



K 



r(Z N {s)f - r(Z N (s-)) 2 - E 2r(Z N ^( S -))(9 i r)(Z JV \s-))AZ? \s) 



i=l 



K K 



= E E(^ r + Tdidjrmbz? ( S )Azf ( S ) 

for some £ such that Z (s— ) — £|| < || AZ iV (s) || < jj, so, as before, there is a constant C such 
that this is bounded above by 

^^^ cEE|AZf(,)||AZf(,)| <CA-^|AZf (.) 

i=l j=l i=l 

Thus, since Z^ is a quadratic pure-jump process, 

K 

e 2NX+(s-t)\ A7 N (c ^ 

0<s<t i=l 



'A^Al 2 



l^wi <c^ e E e2ArA+(s " t) l AZ f( s )r = c '^Er e2iVA+(s ~ t)d [ M ^] 

0<s<t i=l t=l ^ 0+ 



IS 



and, as above, 



for some constant C . 



E 



sup \e N (t)\ 

t<TAT S 



<C'[ I e 2JV * + (-*> ds ] . 



Combining all the above, we see that 



E 



SUp T 2 (Z N (t)) 
t<TAT S 



°u 



9.5. Proof of Theorem [2J Theorem [2] is a consequence of Theorem 5.4 in [TJ; we begin Ito's 
formula to 7Tj(Z Af (t)), from which we obtain an SDE for II (£). 

(29) 7T l (Z N (t))=7T t (Z N (0))+Y] / (0^0(7%-)) dZf (a-) 

i=l ■ / °+ 

+ 5 E f t (P i a j n)(Z lf {s-))d[Zi T 1 Zf] a + r^{t) 
l<i,3<K + 

= TTiCZ^CO)) + JV /"* (Dtt,)(Z 7V ( S -)) • F 7V (Z iV ( S -)) ds + f] f (d m )(Z N (s-)) dM^s) 

K ,t 
+ £/ (^7T i )(Z iV ( S -)MMf] s + £ JV (t) 



where 



(30) rf(t) = £ {7r J (Z^( S ))-vr,(Z^( S -))-£^(Z JV ( S -))AZf( S ) 

0<s<i ^ i=l 

~\ E (^9^ J )(Z JV ( S -))AZi v ( S )AZf( s )}. 



i<«,j<^ 



We consider each of the components individually: 



9.5.1. Finite variation terms. The first three lines in (jlip arise from considering those terms in 
(|29p arising from the finite variation component, 



N [ (Din)(Z N (s-))-F N (Z N (s-))ds. 
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By definition, we have 



NF^{Z lv (s-)) = N 



£ n 4 /3f n (Z%-))Zf( S -)-^(x)Zf( S - 



neNjf 



Gathering terms appropriately, this equals 



K 



i=i 



3=1 



(32) 
(33) 



+N 



/3f (Z%-)) - MZ N {*-))) - (5f(Z N (s-)) - *(Z*(,-))) 

+NF i (Z N (s-)). 



Now, recall (|19p . we have 



N(z"H) 



A(7T(Z^( S -))) 



DF)(tt(Z jv (s-))) e" J o+ 



-/ o 7 (zJV( - )) (Dff)(0.Z^(5-))«fa 



J> 



Theorem [T] shows that r(Z (s— )) — )■ 0; taking the limit as A/" — > oo on both sides gives 



(34) 



(D7T)(7T(Z( S -))) 



1 



■(DF)(7r(Z( fl -))) , 



A(tt(Z( S -))) 
where we have used Corollary [1] to obtain convergence of Z N (s— ) to 7r(Z(s— )). 



First, by assumption, Nfx^Ax.) — > 6ij(x), which combined with some simple manipulations show 
that multiplying (|3ip by (diiTi)(Z N (s— )) and summing over i yield the first two terms in (jlip . 



In a similar fashion, observing that (J32J) converges to 

a(7v(Z(s-))), 
we obtain the third line in the generator pip . 



It remains to show that (|33p vanishes. To that end, we observe that 

(Dvr.XZ^-)) . F N (Z N (s-)) = e, • ((D 7 r)(Z 7V ( S -))F 7V (Z Af ( S -)) , 
while from (|19p . we have 

1 



N/rjN, 



(D7r)(Z JV (s-))F JV (Z JV (s-) 



A(w(Z"(*-))) 



CDF)( W (Z" (*-))) 



e -/-n^ ls - JJ (^(^^ (s -))d Sgr(z Ar (s _ ))i?(z Ar (s _ )) _ 



Now, d22j) yields 



e - /o7 - — (d^z^-))^^^ = y( _ T(zJV(s _ ))Z ^( s -)) = <F0r(Z"(.-))) 
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Thus, recalling from ([8]) that A(x) = (Di?)(x) • ^(x) we have 



1 



A(*r(Z"(a-))) 
r 



(DP)( ff (Z fl H)) )7{n(Z{8-))) 



N, 



A(ir(Z"(*-))) 



y(7r(Z^( S -)))(Di?)(7r(Z^( S -))) <F(tt(Z( S -))) 



^>-)))- "^y -a 



9.5.2. Quadratic variation terms. We next consider the terms 

£ /' (^7rO(Z^(«-))d[J^] a . 
We begin by recalling from (??) that 



[M t N 



* J* AT2 






while rescaling and rearranging ([T]) gives 

1 

AT2 



*&(M) = l^n(Z N (s))Z?(s)ds + 1 ^B? n (Nt). 



0+ 



Moreover, 



E 



_1 R^ 



B&(Nt) 



N 4 



E 



B 



N 



Nt 



^[ B HM = 4 f ^(z^WJz/^Wcfa, 



iV 4 



AT2 



0+ 



so 



^pB&iNt) - J Q+ ^ n (Z N (s))Zi f (s) ds -+ 
in probability as N -> oo. Similarly, ^Df(Nt) - f* + 6^(Z N (s))Z? (s) (is -»• 0. 

Combining these with assumptions ([3]) and @, it follows that as iV — > oo, 



(35) 



H N ]t^ [ A(*-(Z(s))+^(7T(Z( S )))d S . 
J0+ 



Lastly, noting that for x € fi, <5j(x) = /3j(x), we have 

P- lim [Mf] t = f h{<Z{s))+Mir{Z{s)))ds. 

N^oo 1 J Q+ 



Lengthy, yet essentially rote, calculations using (??) and f|35[) show that in the limit as N — > oo, 



J2 (d?ir t )(Z N (s-))d[M q 



fl. 



gives the fourth through eighth lines in (lllj) . 
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9.5.3. Martingale terms. Then, 



M 



N 1 
n,b\ 



t d 



B N 1 

i.nei 



Ns 



N 2 J 0+ W ne ,(Z N (s))Z»(s) 



1 /" dBtl(Ns 



which, using (pQ), equals 



N 2 Jo + ^ nei (Z N (s))Z?( S y 

dB^XNs) 



Jo+ iV Jo 



o+ ' Nj, + ^(ZN {s))Z N {s y 



Now, 



E 



1 f l dB"(N 8 



Nj 0+ P" ne (Z"(s))Z»( S ) 



N 2 



E 



d 



B 



N 



iV« 



° + [^,n e X^{s))Zf{s) 



1e 

iV 



0+ /3f ne ,(Z^( S ))Zf( S ) 



so 



J, 71,0 



■>• £ as N — > oo. Since W^ b (t) is a martingale with jumps tending to as N — > oo, 

it follows from the Martingale central limit theorem that W^ b — > Wi tU: b, where the W^ n & are 
independent Brownian motions. Similarly, the processes W^ d converge in distribution to Brownian 
motions W^a. 



Rearranging, we observe that 



1 

N 



?A ; 



BiU(Nt) 



and 



^Df(Nt) 



o+ 



\l^ nei ^{s))Zf{s)dW^ h {s) 
^5»(Z»(s))Z t N ( s )dW» d (s). 



Lastly, for n ^ |n| ej 



N 



Bi, n (Nt) 



E 



N/ 



fiL(Z»(8))Z?(8)d8 



which, by Assumption (jSJ) is uniformly O(^) on compact sets. Thus, 



jjB itn (Nt) A 0. 



9.5.4. Remainder terms. Lastly, from Taylor's Theorem, we have 

K 



1 f 1 

7r fc (x + Ax)-7r fc (x) -^2(diTr k )(x.)Axi(s) = - ^J Ax t Axj (l-t) 2 (d i d j TT k )(yi + tAx 



eft 



i<*,i<^ 
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so that 



W= E \\ E AZf( S )AZf( S ) / (l-t) 2 (^ i 7r fc )(Z 7V ( S -)+iAZ iV ( S ))-(^7r A; )(Z Ar ( S -))(it 



77 

0<s<t v " l<i,j<K 



(*)| < E (? E AZf(,)AZf( S ) sup |(a^ fc )(Z^( S -)+h)-(^vr fc )(Z w ( S -))|). 

n^„<* ^ Z 1<-,' ;<-&- {h:||h||<i} J 



and 

1^ . _ _ 

0<s<4 *- Z l<i,j<K {h:||h||<-i} 

From our expression (??) for (c?jdj7Tfc)(x), it can be seen that this function is locally Lipschitz 
continuous, so 

\n N (t)\«^[Z% = ^[M%, 
so 77 (t) weakly converges to as N — > 00, by arguments similar to those above. 



9.5.5. Well-posedness for the Martingale problem. 



9.6. Proof of Corollary [21 Let en = l B i ^r; where 0<r'<r<s<l and B are as in 
Propositions Q] and [2j Then, the aforementioned propositions tells us that 

P{sup \Z N (t)-^ Nt Z N (0)\ > N- r } < N~ r , 

t<EN 

while 

F{t(Z n (e n ) > N~ r '} < N~ r . 

In particular Z (e/v) — ^ 7r (Z(0)) € 0. Now, consider the process 

Z Af (t) = Z 7V (t + e^). 

By the above, 2^(0) -^ Z(0) G O, while Z^ satisfies (J23D - Thus, we may apply Theorem |2J to 
conclude that Z^ — >■ Z, where the limiting process has generator (jlip . Moreover, 

Z"(i + £jv ) - Z"(t) = 4> Nt+N£N Z N (0) - ir(Z N (0)). 

Since -/Vejv -> 00 as iV -> 00, the latter converges to uniformly in t, so we must have Z N (. + 

£N ) -z^Ao. 

Lastly, define 7^ : R + ->• M+ by 

iV^-v def f 0<ift<£Af 

I t — £iv otherwise. 



Then, 7 (t) — > £ uniformly in TV, while 



Ar . N ... Jz^o) 0<ift< ejv 



Z iV (7 iV (t))- ,* N 



Z (£) otherwise. 
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Now, the map (/, g) — > f o g ; Dq[0, oo) x Djj+ [0, oo) — > Bn[0, oo) is continuous at points where the 
pair /, g are continuous. Since Z is continuous, we may apply the continuous mapping theorem to 
conclude that Zi N o ^y N — > Z as N — > oo; on the other hand, 



sup 
t 



7»{t)-Z«(>f(t)) =sup Z*(t)-Z N (e N : 

*<£JV 



which converges uniformly to 0. Thus Z — Z — > 0, and the result follows. 
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